Ab initio calculation of H + He+ charge transfer cross sections for plasma physics 
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The charge transfer in low energy (0.25 to 150 eV/amu) H(nZ) + He^(ls) coUisions is investigated 
using a quasi-molecular approach for the n = 2, 3 as well as the first two n = 4 singlet states. 
The diabatic potential energy curves of the IfelT'" molecular ion are obtained from the adiabatic 
potential energy curves and the non-adiabatic radial coupling matrix elements using a two-by-two 
diabatization method, and a time-dependent wave-packet approach is used to calculate the state- 
to-state cross sections. We find a strong dependence of the charge transfer cross section in the 
principal and orbital quantum numbers n and I of the initial or final state. We estimate the effect 
of the non-adiabatic rotational couplings, which is found to be important even at energies below 
1 eV/amu. However, the effect is small on the total cross sections at energies below 10 eV/amu. 
We observe that to calculate charge transfer cross sections in a n manifold, it is only necessary to 
include states with n < n, and we discuss the limitations of our approach as the number of states 
increases. 

PACS numbers: 34.70. ~(-c, 52.20.Hv 



1. INTRODUCTION 

Starting with the historical work of Massey and Smith in 1933 [l| and due to the apparent simplicity of this two 
electron system, the asymmetrical charge transfer process He (Is^ ^5') -I- H"*" — > H -I- He"'" has quickly been considered 
as a prototype for semi-classical methods to treat collisions 0-0 • This process is dominated by the capture into the 
H(ls) state but later on, the charge transfer excitation and the direct excitation processes have been studied in detail 
both theoretically and experimentally ,7.-.15.j. All these reactions require intermediate collision energy which suits 
perfectly the semi-empirical description of the collision. It is also the case for the charge transfer mechanism that 
involves H" ion as projectile (He^""" -I- H" -J> H -I- Y{e+{nl)) ^.\vf\. 

At much lower collisional energy, charge transfer can populate excited states of He(lsn^ ^''^L) from the corresponding 
excited states nl of H. Although there is no measurement for those processes, they have been studied theoretically 
using a semi-classical approach with a linear trajectory for the nuclei taking into account the coupling between the 
final Stark splitting H states and the initial He states at large internuclear distances {R > 20 a.u.) and neglecting 
electron translational factors and rotational couplings (isl . [l9j . These works provide data for an energy range between 
2.5 eV/amu and 10 keV/amu. 

Recently, it has appeared that those low energy charge transfer processes involving excited hydrogen states could 
be of major importance for the monitoring of warm plasmas [20| . Indeed, spectroscopic methods are among the 
most effective approaches to determine particle transport in magnetic confined fusion plasmas. From this point of 
view, simulations of excited He emissions resolved in space and time has been proposed as a tool independent of the 
theoretical plasma model. Atomic physics simulations can be compared to experimental data and the diffusion and 
convective velocity parameters could be determined. However, the He ions will interact with the H/D background 
via a charge transfer mechanism which modifies the population of the He excited states and therefore the intensity of 
the emission lines. A self-consistent approach to the description of the coupling of the radiating He with the plasma 
background via charge transfer has shown that these processes are very important at low collisional energies (typically 
of the order of 0.1 eV to 100 cV) and that, therefore, an accurate knowledge of charge transfer cross sections in this 
energy range is essential [20| . 

In this work, we have adopted a quasi-molecular approach of the ion-atom collision based on the use of quantum- 
chemistry ab initio methods to obtain the potential energy curves as well as the radial and rotational coupling matrix 
elements of the quasi-molecule HeH"*". A wave packet method is used to treat the curve-crossing dynamics resulting 
from the failure of the Born-Oppenheimer approximation. A Gaussian wave packet is prepared in the entrance channel 
and propagated on the coupled ro-electronic channels. The collision matrix elements are computed from an analysis 
of the flux in the asymptotic region by using properties of absorbing potentials, giving access to the charge transfer 
cross-sections for the processes He"*" (Is) -I- H(n^) He {Isn'l' ^-^L) + H+ where rt = 2 — 3. We estimate the influence 
of the rotational couplings on the cross section. 
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2. THEORY 
2.1. Molecular data 

The Hamiltonian is given as the sum of an electronic part and a nuclear kinetic part: 

H = T^ + H''^ (1) 

The electronic Hamiltonian includes a kinetic term for the electrons and all the potential energy terms. 

The potential energy curves (PEC) UmA and the adiabatic electronic functions (^„iA solve the electronic motion: 

H'^'CrnAir; R) = t/™A(i?)CmA(r; R) , (2) 

where r stands for the electron coordinates and R is the radial coordinate for the nuclei, m is used to number the 
states for a given A, which is the quantum number associated to Lz, the projection of the total electronic orbital 
angular momentum L onto the molecular z axis. The molecular electronic states are classified according to the value 
of |A|: S states correspond to A = 0, 11 states to |A| = 1, and A states to |A| = 2. We therefore see that states with 
|A| 7^ are doubly degenerate for singlet states. In the atomic limit (R — >■ oo), A becomes mi, the magnetic quantum 
number. 

On the other hand, can be written in atomic units as the sum of a radial part, 



1 

2^ 

where fi is the reduced mass of the system, and a rotational part given by 



i/rad ^ __g2 ^ (3) 



2fiR^ 2fiR^ 



K^ + L^ -2KzLz-K+L^-K^L+ (4) 



where N is the nuclear angular momentum. In this work, we will focus on singlet states and we will not consider 
any spin-dependent interactions, so that the total angular momentum is K = N + L. Since = jjiad _^ jj'o^^ the 
nuclear wave function is the product of a radial part and an angular part: ipmACR) — '^PmAiR) \KAM). The angular 
functions are eigenfunctions of and with eigenvalues K{K + 1) and A, respectively. The action of the ladder 

operators K± = K^± iKy is given by K±\KKM) = [K{K + 1) - A(A =F 1)] ^''^\KK T IM). 
In the basis of these electronic-rotational functions, the matrix elements of i/'"* are given by 



H'°lK,m'A'K' = ^1 {k{K + 1) -K^)5,nm'5AA' - {L-Uw [K{K + 1) - A(A - 1)]^/' 5a,A'-i 



-(L+)„,„,- [K{K + 1) - A(A + 1)] ^'^5aa'+i \5kk' (5) 



where the contribution from (L^ + L^)mm'^AA' has been neglected. We see that states with AA — ±1 will interact 
through the rotational Hamiltonian. 

To treat the effects of the rotational Hamiltonian, it is more convenient to work with parity adapted functions [2H . 
These functions are defined by 



Imi^AMe) ^ 



v2 + 2Sao 



|i^AM)C™A + {-ir(-\K - AM)Cm-A (6) 



where e = 1 and e — —1 correspond to e and / states, respectively. Using ([5]) and ([5]), it can be shown that only 
connects states of the same parity. As states can have e or / symmetry, only half the H states must be taken 
into account in the calculations. 

We will consider here the n = 1 — 3 ^S"*", ^H and ^A states as well as the first two n = A states. We 
could not include more n = 4 states in the calculations, as we were not able to calculate the radial non-adiabatic 
couplings between these states by ab initio methods. The dissociative atomic states and the asymptotic energies of the 
molecular states are shown in table H] We have also considered the n — 2 and '^H states to allow the comparison 
with [3, [1^ (see section [3761) . The adiabatic potential energy curves (PEC) for these states have been calculated 
using the ab initio quantum chemistry package MOLPRO version 2006.1 [l^]. An adapted basis set consisting of the 
aug-cc-pv5Z basis set psj supplemented by one contracted Gaussian function per orbital per atom up to n = 4 has 
been used. Details of the calculations can be found in [23|. Due to the electric field produced by He+(ls), there 
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TABLE I: CASSCF energies at i? = 70 a.u. and dissociative products of tlie singlet states included in the calculations. 
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is a Stark mixing of the hydrogen states. However, these Stark states adiabatically become pure atomic states as 
i? — > oo so that we know with certainty the atomic configuration corresponding to a molecular state, as indicated in 
table m The PEC have been calculated at the state-averaged complete active space self-consistent field (CASSCF) 
level and are shown in figure [T] for the n = 2 — 3 states. This approach allows us to compute the radial non-adiabatic 
coupling matrix elements -F'mA,m'A — (CmAl^fllCm'A)j which are used to build the diabatic representation [25j. The 
adiabatic-to-diabatic transformation matrix D is the solution to the differential matrix equation ds^ -I- F • D — 
and the diabatic potential energy curves are the diagonal elements of the matrix U'^ = D"-'^ • U • D, where U is the 
matrix of H'^^ in the adiabatic representation. We have used an approximate F matrix by keeping only the couplings 
between adjacent states, i.e. the elements Frn,m+\. The main reason is that this approach simplifies considerably the 
diabatization procedure [23|. It has been shown to give similar results as when the complete F matrix is taken into 
account in dynamical calculations (26j , something we have also observed in the low-energy calculation of the charge 
transfer cross sections for the n = 2 states (see below). 

The non-adiabatic rotational coupling matrix elements (CmAl-^-i-ICm'A/) appearing in equation ([S]) have also been 
computed at the CASSCF level using MOLPRO. As pointed out in [23,|2g| and 2J|, some of these couplings behave 
asymptotically as i?, which is due to the choice of the electronic coordinates. Due to the factor XjB? in equation ([5]), 
these couplings will thus decrease as 1 /i?, much slower that the radial couplings which decrease to zero extremely fast 
outside the interaction region. This causes a problem in the calculation of the cross sections as it implies the use of 
very large numerical grids that increase tremendously the calculation time. To solve this problem, we modified the 
problematic rotational couplings outside the interaction region where we required that they decrease to the atomic 
values of the couplings. We have tried various switching functions to find a set of parameters that had no effect on 
the cross sections. This approximation is also justified in our case by the fact that the linear rotational couplings 
usually connect two states in the same atomic configuration, so that the modification will not influence the charge 
transfer cross sections. 

There are also cases where the atomic value of (CmA|i±|Cm'A') is a constant, but not zero. The rotational Hamil- 
tonian then decreases as which still implies the use of large numerical grids. However, this can only happen 

again for transitions between two states in the same atomic configuration (electron excitation), a process we do not 
consider here. 
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In the atomic limit, when the non-adiabatic rotational couplings are neglected, our method implies conservation of 
the magnetic quantum number m^^ as the Hamiltonian is diagonal in A. When they are included in the calculations, 
we have interaction between states with Am^ = 0, ±1. 
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FIG. 1; Adiabatic (left) and diabatic (right) PEC of the n = 2, 3 states of HeH+. Full lines, states. Dots, ^11 states. 
Dashed lines, states. The n = 1 ^'E^ states have been excluded from the diabatization procedure. 



2.2. Cross section calculation 



The cross section corresponding to the transfer of an electron from an initial state m, A to a final state to', A' is 
given by 

<Jm'A',^A{E) = -^^Y.(^K + 1)|5,^,a,,„,a(£;) - 5,n'mSA'A\^ (7) 
'^7nA\-'^) K 

where k„iA is the wave number in the entrance channel, kmA — •\/2/i(£' — UmA)- As the Hamiltonian is diagonal in 
K (see equation dS])), the cross section must be calculated for each value of K until convergence in equation ([7]). 

We use the coupled channel formalism in the rotational-electronic diabatic representation. In the time-dependent 
formalism, we start by defining a Gaussian initial wave packet which is propagated in time using the split operator 
algorithm |30]- The coupled equations give access to the wave packets on all the rotational-electronic states. For each 
value of K, the scattering matrix elements I'S'^'A'.mAl^)!^ then extracted using the flux operator formalism with 
a complex absorbing potential f3l|, [s^ . 

We start by defining the functions 

^^A,E = y^^4(fc™A^)C™A (8) 

where h^{kmAR) are the Riccati-Hankel functions [33j and CmA ^he electronic wave functions in the diabatic 
representation. 
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We then introduce the time- independent energy normahzed wave functions I^'^'a^b); sohitions of 

JoA|,T, + ,A' \_TT . I.Tr + .-'f 



Yy^^" + VR^^^(^^ + " ^ ^ + ^ ) ^ u^^\'^t;lE) (9) 

and satisfying the asymptotic condition 

, A' 

These stationary eigenfunctions can be constructed as the Fourier-transform of a time-dependent wave packet 



l*mf - TT^TTTn / l*™A(i)) exp(zM) dt . (11) 



The vector is constructed by propagating an initial wave packet |^(0)) in time using the Hamiltonian matrix 

in the rotational-electronic diabatic representation: 

=exp(-ffl^t)|*(0)) . (12) 

The initial wave packet is zero except in the diabatic channel mA, where it is represented by a Gaussian function 
g{R) of width a and centered around Rq: 



1 



liR) = I exp 



ikoR 



{R - Ro) 



(13) 



rjj^ is the amplitude of the initial wave packet on the stationary states: 



r^A = (KuaJ^o) = h+{k^AR)9{R)dR . (14) 

The flux operator is defined by [3l| 



where Rc is a point in the asymptotic region {i.e. such that there is no interactions for R> Rc) located behind ido- 
lising equations pUj) and (ITSt . one arrives at 

{"^t^lEWK^E) = ^ E \S^^'^',mA{E)? (16) 
m',A' 

The sum can be removed using the projector onto the electronic state m'A', Pm'K' — ICm'A')(CTO'A' I' obtain 

{<:lE\P^'^'FP^n'A'\^t;lE) = ^l^^'A',™A(i?)P (17) 

A complex absorbing potential ~iW is then added to the Hamiltonian which becomes H' — H ~ iW. The left hand 
side of equation (|17p is then calculated using equations (ITT|) and ([T^ with H' instead of H. This is allowed if the 
CAP is "switched on" in the asymptotic region: in the interaction region, the CAP vanishes and the value of $(i) 
propagated with H or H' will be identical. Combining equations (ITT|) and pT)). we find that the state-to-state cross 
section is given by 



oo poo 



\S^'A',^AiE)\^^ ^ ,m\2 / dt'{m\Pm-A'WP^,K,m'))eMm' ~t)] (is) 

^t^W^aKE)V h Jq 

Equation (flS]) is used to obtain the matrix elements of S. 
In our calculations, we used a CAP given by 

mR) = vc '^^~^f (19) 

where rjc is the strength of the CAP and i?oo is the last point of the grid. 
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3. CHARGE TRANSFER CROSS SECTIONS 
3.1. Computational details 

3.1.1. Parameters for dynamics 

The parameters are chosen so as to ensure convergence of the sum in equation ([7]) while keeping the norm of the S 
matrix close to unity. For the calculations involving the n = 2 states, a typical set of parameters consists of 2^^ points 
for a grid of 60 a.u., an initial wave packet located around i?o — 40 a.u. of width a — 0.2 and a CAP starting at 
Rc = 45 a.u. of strength rjc — 0.01. The time needed for the wave packet to return to the asymptotic region obviously 
depends on the collision energy. For the set of parameters above, it is approximatively contained between 2 • 10^ and 
3 • 10^ a.u. for energies between 100 and 0.2 eV/amu. 

For the calculations of the cross sections involving n = 3 states, we had to use grids up to i? = 100 a.u. This is 
due to the fact that the number of avoided crossings increases strongly with n, so that the positions of the radial 
non-adiabatic couplings are shifted to larger internuclear distances [2J|. The time needed for the wave packet to 
return to the asymptotic region is therefore increased, and can be as high as 5 • 10^ a.u. for low energies. 

When the rotational couplings are included in the calculations, the convergence of the partial cross sections is 
considerably slower as a function oi K . It is therefore necessary to use much larger grids and the time of propagation 
is therefore increased tremendously. In addition, the number of points of the grid must also be increased to keep a 
constant step dR, again extending the computational time. In order to reduce the calculation time, we used grids of 
variable size ranging from 150 a.u. for small K to 600 a.u. for values of K around 2500. When treating the n — i 
states, we could not perform calculations of partial cross sections for energies higher than 10 eV/amu as even these 
large grids did not ensure convergence. The width of the wave packet was also increased up to cr = 1.5 for these large 
grids, as wave packets with larger width stay more compact. 

3.1.2. Non-adiabatic radial couplings 

The first result that needs to be established is the validity of our approximation which consists in only retaining the 
non-adiabatic radial couplings -Fm,m+i instead of the complete F matrix. We show in figure [2] a comparison between 
the two methods in the calculation of the cross sections for the process H(2s) + IIe+(ls) H+ + He(ls2Z ^L). We 
conclude that this approximation is perfectly valid at low energies but that small deviations are observed at higher 
energies, E ~ 100 eV/amu. The same conclusion is reached for the cross section for the process ii{2p) + He"'' (Is) — > 
H+ + He(ls2; iL). 

25 
20 
15 
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FIG. 2: Comparison of the two-by-two diabatization with the use of the complete F matrix in the calculation of charge transfer 
cross sections with He'*"(ls) -I- H(2s) in the ^E"*" symmetry as the initial state. Full line and -|- signs; charge transfer onto 
He(ls2p ^P°) + with two-by-two or complete F matrix, respectively. Dashed line and x signs: charge transfer onto 
He(ls2s ^S) + with two-by-two or complete F matrix, respectively. 
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Another issue is the fact that due to the Stark effect on hydrogen, some PEC undergo avoided crossings at large 
internuclear distances. However, as was pointed out in [23|, the large amplitude and the narrowness of the non- 
adiabatic radial couplings at those points indicate that a full diagonal diabatic representation at the crossing is 
perfectly justified, so that these crossings will not affect the cross sections. 

Finally, as pointed out in 24] , the effect of the electron translation factors was found to be negligible in the range 
of energy considered in this work. 



3.2. General observations 



From a practical viewpoint, the calculation time of a cross section goes roughly as e*^ "^™, where m is the number of 
states included in the calculations, so that the computing time doubles every time two additional states are considered. 
This again extends the computational time when rotational interactions are taken into account, as E and 11 states 
must be considered in the same calculation. It is therefore important to take the least states possible, and we will see 
that states with different values of n can be considered independently. 

A few things seem to come out from our cross section calculations, which are presented below. The first is that, to 
a good degree of precision, the charge transfer cross sections in a given n manifold are not modified by the inclusion 
of states with a principal quantum number n' ^ n, a fact illustrated in figure [3] for the n = 2 ^E"*" states. It was 
anticipated that the n — I states did not play any role in the n > 1 cross sections since they are much lower in 
energy, but the possibility to treat the n = 2 and n = 3 states separately was much less obvious. Indeed, the shape of 
the diabatic PEC is strongly influenced by the inclusion of states with different principal quantum number and the 
density of states increases with n [13] . We reached the same conclusion for the influence of the first two n = 4 states 
on the n = 3 manifold and in the ^11 symmetry. 

However, while the cross section from the n = 2 onto the n = 1 states is always negligible, this is not the case in 
general as the states interact through non-adiabatic radial couplings. For example, the cross sections from the n = 3 
onto the n = 2 states is not negligible, although the cross sections inside the n = 3 manifold are not modified by the 
inclusion of the n = 2 states. This means that the cross sections onto the n — 2 states result from a decrease of the 
elastic cross section. 
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FIG. 3: Illustration of the possibility to treat states in different n manifolds independently in the calculation of charge transfer 
cross sections with He^(ls) + H(2s) in the ^E^ symmetry as the initial state. Full line: charge transfer onto He(ls2p ^ P°) + 
using the four n = 2 states. + signs: the same, but with the two n = 1 states included. □ signs: the same, but with the 
six n — 3 states included. Dashed line: charge transfer onto He(ls2s ^S) + H"*" using the four n = 2 states, x signs: the same, 
but with the two n = 1 states included. Q signs: the same, but with the six n — 3 states included. 



The second observation is that there is a dependence of the cross section with n. This dependence was expected 
since the cross section scales classically as for Rydberg states. 

Finally, we also observed that, in a n manifold, the cross section always increases with the orbital quantum number 
I of the initial state. 
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3.3. n — 2 states 

The cross sections with H(2p) and H(2s) in the E symmetry as weh as H(2p) in the 11 symmetry are presented in 
figures S) 

The behavior of the cross section for the process H(2p) + He"*" (Is) — ?► He(ls2p) + H"*" is completely different in the 
E (figure [4)d) and 11 (figure [4j;) symmetry: the total cross section from H(2p) will be governed by the 11 states at low 
energy and by the E states at high energy [E > 100 eV/amu). Another difference is the behaviour of the cross section 
when the rotational couplings are included in the calculations: they have no effect for the transition with H(2p) + 
He"*" (Is) in the E symmetry as the initial state, but strongly modifiy the cross section for the corresponding 11 state. 
The cross section between the two 11 states is decreased while the cross section from 11 to E states is increased so that 
the total cross section with H(2p) + He+(ls) as the initial state is roughly the same as when the rotational couplings 
were neglected. 

The transition H(2s) + He"*" (Is) — )• He(ls2p) + H"*" is also affected by the inclusion of the rotational couplings at 
energies i? > 10 eV/amu (figure|4^). However, in this case the cross section between E states and from E to H states 
are both increased. This simply means that for this state, the elastic cross section is decreased by the inclusion of 
rotational couplings. 

3.4. n = 3 states 

There are 12 n = 3 states. The cross section between the six -'^E"'" states, presented in figures [SI H] and [71 have been 
calculated including the n = 2 states since the cross section from the n = 3 to the n — 2 states is not negligible, as 
shown on these figures. There are also four 7i = 3 *^H states and two states. The charge transfer cross sections 
between these states are not presented here, but are available upon request to one of the authors, along with all the 
cross sections presented throughout this article. From these figures, it is clear that there is a dependence of the cross 
section in the principal and orbital quantum numbers, n and I, of the initial H(nZ) state: the charge transfer cross 
section is much larger for n = 3 than for n = 2, and also increases with the value of I. The other difference between 
the n — 2 and n = 3 manifolds is the influence of rotational couplings. While in the n = 2 states they were not 
influential at energies _E < 10 eV/amu, this is not the case for the n = 3 states where they play an important part 
even at energies below 1 eV/amu. We observe the intuitive fact that the cross sections between E states are smaller 
when the rotational interactions are taken into account, corresponding to the fact that a part of the cross section is 
transferred onto the H states. 

We also observe that the cross sections onto the He(ls2p ^P°) + H+ state in the E symmetry is smaller when the 
initial state is higher in energy. The cross section to the other n = 2 states, Hc(ls2s ^S) + H"*" in the E symmetry 
and He(ls2p ^P°) + H"*" in the H symmetry, are negligible and therefore not shown. Interestingly, the charge transfer 
cross sections from the n = 2 states onto the n = 3 states are all negligible. 

3.5. n = 4 states 

We have included the first two n = 4 singlet states in the E symmetry. We could not consider more than the 
first two n — 4: states in the diabatization since we were not able to calculate the radial non-adiabatic couplings for 
the higher-lying states. We calculated the cross section starting from the second n = 4 state, H(4p) + He"'"(ls). It 
confirms once again the dependence of the cross section in the quantum number n for a given value of I. We also 
observe that the cross sections from this n — 4 state onto the n = 3 states are not negligible, as shown in figure [S] 

Moreover, it is observed (not shown) that despite the fact that the n = 3 and n = 4 states are close in energy and 
interact through radial non-adiabatic couplings, the cross sections in the n — 3 manifold are not infiuenced by these 
two 71 = 4 states, with the exception of the cross section with H(3s) + Hc+(ls) as the initial state, which is slightly 
modified at energies higher than 10 eV/amu. In addition, the cross section from H(3s) into the first n = 4 state, 
He(ls4s) + H+, is negligible. We thus reach the same conclusion as in the case of the n — 3 states: the cross sections 
from the n = 3 states onto the n = 2 states are not negligible, although they do not influence the cross section within 
the n = 2 manifold. 

In conclusion, it would thus seems that the calculations of the charge transfer cross sections within a given n 
manifold require to take into account only the states with n' < n. 
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Energy (eV/amu) 

FIG. 4: Charge transfer cross sections between the n = 2 states. 

(a) With He+(ls) + H(2s), ^E+, as the initial state. Full hnc: charge transfer onto He(ls2s ^5") + H+, ^S+. + signs: same, 
but with rotational couplings. Dashed line: charge transfer onto He(ls2p ^P°) + H"*", ^E^. x signs: same, but with rotational 
couplings. □ signs: charge transfer onto He(ls2p ^P°), ^11. 

(b) The same as (a), but with He"'"(ls) + H(2p), ^'S'^ , as the initial state. 

(c) With He+(ls) + H(2s), ^n, as the initial state. Full line: charge transfer onto He(ls2p ^P°) + H+, ^U. + signs: same, but 
with rotational couplings. □ signs: charge transfer onto He(ls2s ^S) + H"*", ^E"*". Q signs: charge transfer onto He(ls2p ^P°) 
+ H+, iS+. 
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Energy (eV/amu) 

FIG. 5: Charge transfer cross sections with He^(ls) + H(3d), ^E^, as the initial state. Full line: charge transfer onto 
He(ls3s ^S) + H^, ^E"*". + signs: same, but with rotational couplings. Dashed line: charge transfer onto He(ls3d ^D) + H^, 
^E+. X signs: same, but with rotational couplings. Light dashed line: charge transfer onto He(ls3p ^P°) + H+, ^E"'". □ signs: 
same, but with rotational couplings. Dots: charge transfer onto He(ls2p ^P°) + H+, ^E+. Q signs: same, but with rotational 
couplings. A signs: charge transfer onto He(ls3d ^ D) + H'*', ^11. V signs: charge transfer onto He(ls3p ^ P°) + H^, ^11. 




FIG. 6: Charge transfer cross sections with He"'"(ls) + H(3p), ^E^, as the initial state. Full line: charge transfer onto 
He(ls3s ^S) + H^, ^E""". + signs: same, but with rotational couplings. Dashed line: charge transfer onto He(ls3d ^D) + H^, 
^E^. X signs: same, but with rotational couplings. Light dashed line: charge transfer onto He(ls3p ^P°) + H+, ^E"*". □ signs: 
same, but with rotational couplings. Dots: charge transfer onto He(ls2p ^ P°) + H"*", ^E^. Q signs: same, but with rotational 
couplings. A signs: charge transfer onto He(ls3ci ^D) + H^, ^11. V signs: charge transfer onto He(ls3p ^ P°) + H+, ^11. 



3.6. Total cross sections 



The total cross sections starting from a given nl state of H are obtained by summing all the contributions from 
within a A manifold 

oo n' — 1 

a{nlK) = XI 51 '^("'^ ^ "'^'^) (20) 

n' = l i'=0 

and then by summing the contributions from all the A [l9j : 

A=-i 

It should be noticed that as the states with A 7^ are doubly degenerate, they contribute twice in the sum in equation 
([21]) . When the rotational couplings are taken into account, the only difference is an additional sum over A' in equation 

(uni)- 
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FIG. 7: Charge transfer cross sections with He^(ls) + H(3s), ^E^, as the initial state. Full line: charge transfer onto 
He(ls3s ^S) + H+, 1E+. + signs: same, but with rotational couplings. Dashed line: charge transfer onto He(ls3d ^D) + H+, 
^E+. X signs: same, but with rotational couplings. Light dashed line: charge transfer onto He(ls3p ^P°) + H+, ^E"'". □ signs: 
same, but with rotational couplings. Dots: charge transfer onto He(ls2p ^P°) + H+, ^E+. Q signs: same, but with rotational 
couplings. A signs: charge transfer onto He(ls3d ^D) + H^, ^11. V signs: charge transfer onto He(ls3p ^ P°) + H+, ^11. 




FIG. 8: Charge transfer cross sections with He'''(ls) + H(4p), ^E^, as the initial state. Full line: charge transfer onto 
He(ls4s ^S) + H"*", ^E^. + signs: same, but with rotational couplings. Dashed line: charge transfer onto He(ls3p ^P°) + H"^, 
'^E^. Light dashed line: charge transfer onto He(ls3ci ^D) + H^, ^E""". Dots: charge transfer onto He(ls3s ^S) + H"*", ^E^. 



The total cross sections from the H(nZ) states, with n = 2,3 are shown in figure [5] on a log-log scale. This figure 
clearly illustrates the dependence of the charge transfer cross section in n and I . It also shows that while the inclusion 
of rotational couplings modifies the behavior of state-to-state cross sections at low energy, it modifies only slightly 
the total charge transfer cross section of a n manifold. At energies higher than 10 eV/amu, the influence of rotational 
couplings starts to be important and the total cross section is increased. It would therefore be interesting to investigate 
the contributions of the rotational couplings at higher energies for n — 3 states, but it is clear that our method is not 
adapted to such calculations. 

In the same way, we can determine the total cross section with lle{lsnl ^L) + H+ state as the final state. In 
helium-based plasma diagnostic, a correct estimation of the populations of the various He(lsn/ ^^^L) levels, which 
are modified by charge transfer, is necessary. States such as lle{lsnp ^P°) decay radiatively to the ground state and 
these emission lines can be observed. In figure [TUl we have grouped together the charge transfer cross sections with 
He(ls2p i'3po-j _|_ JJ+ ^j^j He(ls3p ^P°) + 11+ as the final states. Interestingly, the conclusions are similar when 
we consider the sum of all cross sections into a specific final state as when we considered a specific initial state: we 
again see a dependence in n, which is now the principal quantum number of the helium atom. We also see that the 
influence of the rotational couplings is weak at low energies but that the total charge transfer cross section is increased 
at energies > 10 eV/amu. 

These results can be compared with those of Chibisov et al [l^. To describe the charge transfer process, these 
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FIG. 9: Total charge transfer cross sections starting from ¥L{nl) + He'''(ls). 

+ signs: H(2s), x signs: H(2p), Full line: H(n — 2), and ■ signs: H(n = 2) with rotational couplings. 

O signs: H(3s). A signs: H(3p). signs: H(3d). Dashed line: H(n — 3) and A signs: H(n = 2) with rotational couplings. 



authors used a semi-classical method with an atomic basis where only the Stark couplings between the atomic states 
are taken into account and without rotational couplings. The cross section with He(ls2p) as the final state ([l^, fig. 
3) can be compared in the range of energy between 2.5 and 200 eV/amu (see figure [TU]). For the singlet states, we 
see that not only is the behavior at small energies different, but also that the cross section is several times smaller 
in our calculations. In the triplet symmetry, the order of magnitude of the cross sections are roughly the same, but 
the general behavior is different. The comparison can also be made for the n = 3 states for S, 11 and A. In p^ . 
Chibisov et al present state-to-state cross section calculations, so that the comparison with our calculations is direct. 
The results are again qualitatively very different, showing the limitations of their method. 




Energy (eV/amu) 



FIG. 10: Charge transfer cross sections with He(ls2p ^-^p") + and He(ls3p ^P°) + H"*" as the final state and comparison 
with 19]. Full line: He(ls2p ^P°) + H+ as the final state. + signs: same, but with rotational couplings. Dashed lines: 
IIe(ls2p '^P°) + as the final state. Light dashed line: He(ls3p ^P°) + as the final state, x signs: same, but with 
rotational couplings. Circles and diamonds: calculations of jl^ for the singlet and triplet states, respectively. The points were 
extracted graphically. 



4. CONCLUSIONS 



Using a quasi-molecular approach and a wave packet propagation method, we have computed the state-to-state 
cross sections for the charge transfer coUisional process H(nZ) -f He+(ls) — He(lsn'/' ^L') for all the n, n' = 2,3 
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singlet states (as well as the first two n = 4 states in the ^'S'^ symmetry) in the energy range between 0.25 and 150 
eV/amu. We have also investigated the effect of the non-adiabatic rotational couplings on the charge transfer cross 
sections. All the cross sections are not presented in this article, but are available upon request to one of the authors. 

We have found that our method is adapted when the rotational couplings are neglected, but is problematic at 
energies higher than 10 eV/amu when the rotational couplings are included due to very long computational times. 

We have found a strong dependence of the state-to-state charge transfer cross sections in the principal and orbital 
quantum numbers, n and I, of the hydrogen atom. We observed that the rotational couplings have an influence on 
the cross sections even at low energies, but that their effect is increased with n: for n = 2 states, we found that the 
effect of the couplings start to be important at energies higher than about 5 eV/amu, while for the n = 3 states they 
modify the cross sections even at energies below 1 eV/amu. However, the total cross sections are not modified by 
the inclusion of rotational couplings at energies below 10 eV/amu. The effect of these couplings should therefore be 
investigated at intermediate energies by another method. 
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